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ABSTRACT 

In this work, we present the mass and magnetic distributions found in a recent Adaptive Mesh 
Refinement (AMR) MHD simulation of supersonic, super- Alfvenic, self gravitating turbulence. Power 
law tails are found in both mass density and magnetic field probability density functions, with P{p) cx 
p 6 and P{B) cx B~ 2 '. A power law relationship is also found between magnetic field strength and 
density, with B cx p ' 5 , throughout the collapsing gas. The mass distribution of gravitationally bound 
cores is shown to be in excellent agreement with recent observation of prestellar cores. The mass 
to flux distribution of cores is also found to be in excellent agreement with recent Zeeman splitting 
measurements . 

Subject headings: methods: numerical — AMR, MHD 
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1. INTRODUCTION 

Understanding the details of star formation is one of 
the great open problems in astrophysics today. One 
of the central questions is the role of magnetic fields. 
Two opposite paradigms have been explored during the 
study of this problem. The earliest par adigm argues for 
strong fields that do minate the process ([Shu et alj|1987t 
IMouschov"i as 1987ajJ2) and yield their dominance by way 
of ambipolar diffusion; the other extreme argues for rel- 
atively weak fields that at most a lter the thickness o f 
shocks (|Padoan fe Nordlundl [19991 iPadoan et all 12001 , 
with the majority of the details being determined by the 
turbulence. 

Central to the development of the second paradigm 
has been the development of more advance comput- 
ing hardware and software. Only a handful of simu- 
lations have been performed that incorporate all three 
necessary aspects: tu rbulence, gravi t y, and magnetic 
fields. Notably are iGammie et al.1 (I2003T) . iLi et al ' 
(120041) iHeitsch et all (12001ft iVazouez-Semadeni et al 



■2Q^l . lPadoan fc Nordfond ()2009MPrice fc Bate) ([20081 
|2oToF . The enormous range of scales at work in this 
problem have motivated the developme nt of high dy- 
namic range techniques, namely AMR ([Fromang et al.l 
[200l ICollins et al1l20T0h and t he addition of MHD to 
SPH ([Price fc Monaghari 12001 . These techniques are 
beginning to yield intere sting results on the nature of 
MHD in star formation (|Price fc Bate! 120081: iDib et al.l 
I2010I) . 

In this work we focus on the distribution of mass and 
magnetic fields in super- Alfvenic simulation using high 
resolution AMR simulations. Section [2] will introduce 
the numerical algorithm and discuss the initial conditions 
and parameters of the simulation. In Sections [3] and [4] 
we present the volume and mass weighted density PDF 
for both density and magnetic field. In Section [5] we 
discuss the relation between mass density and magnetic 
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field strength. In Section [SJ we discuss the relationship 
between line-of-sight magnetic field strength and column 
density (typically used as a proxy for the mass-to-flux 
ratio M/<3?, where M is the mass, and $ is the magnetic 
flux), and compare to recent observations. In Section [7] 
we present the mass distribution of cores, and compare 
to observations of protostellar cores. In Section[S]we will 
conclude. 

2. NUMERICAL MODEL 
2.1. Simulation software 

The simulation was perform e d with the MHD exten- 
sion of Enzo by ICollins et al.l (12010ft. I t employs the 
MHD Godunov solver of ILi et al.l ([20081 ) to solve the 
MHD eq uations: to constra i n V • B it employs the CT 
scheme of lGardiner fc Stond ([20051 ) ; the AMR re finement 
is don e with the divergence free interpolation of iBalsaral 
(|2001[) . An isothermal equation of state and ideal MHD 
are assumed throughout. 

Refinement is performed such that the Jeans length, 

Aj = 

where c s is the sound speed and G is the gravitational 
constant, is resolved by four grid zones everywhere. 
Thus, whenever the density in a zone is above the "Tru- 
elove Density" 
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that zone is flagged for refinement. The maximum num- 
ber of levels is 4. With a base grid of 128 3 , this gives a 
maximum effective resolution of 2048 3 . 

2.2. Simulation Parameters 

We began with a uniform density and magnetic field 
with no self gravity, and stirred with a Gaussian ran- 
dom field with no compressional modes, with power dis- 
tributed in a top-hat between between 1 < k < 2 and 
mean Mac h number 8.9. Driving is done in the same 
manner as iMac Low! ([1999ft . such that the energy injec- 
tion rate is constant. After a number of dynamical times, 
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gravity was activated. This time is t = 0. Driving con- 
tinued for the duration of the collapse. 

The initial uniform fields have a ratio of gas to mag- 
netic pressure of j3 =8ttc 2 (p) / (B 2 ) = 22.2, but the mean 
squared field strength has been amplified by the driv- 
ing so that when gravity is turned on at t = 0, j3 = 
0.2. This yields an Alfven Mach number of 2.8 at t = 0. 
have been chosen so that a = 5a 2 R/(3GM) sa 1, where 
a is the three dimensional velocity dispersion, M is the 
total mass in the box, and R is half the width of the 
box. The typical parameter used to measure the rela- 
tive effects of magnetic and gravitational energy is the 
mass-to-flux ratio in units of the critical value for sup- 
port. If we use the critical value for a sheet, as in 
INakano fc Nakamural dl978h . (M/$) c = (4tt 2 G)- 1 / 2 , the 
uniform (pre-turbulence) cube gives us 



A = (M/$)/(M/$) c = 18.7. 



(1) 



Mass to flux is a reasonable quantity to parameter- 
ize magnetic support for simple geometries, like spheres 
and disks. However, for situations like the one presented 
here, in both the turbulent box at t = and the prestellar 
cores discussed later, the spatial structure of the objects 
in question are not easily defined. We find it more in- 
structive to examine the ratio of energies. For sake of 
comparison to M/$, we define 



A_e — v Ec/E B , 



where 



E G = G 



E, 



P(r)p(r') 
\r — r'\ 



D 2 



-d 3 r. 



(2) 

(3) 
(4) 



For the initial conditions at t = (after the initial tur- 
bulent evolution, but before gravitational collapse), we 
have 



A K = 6.5 



(5) 



For reference, a sphere with critical mass M c = 
$/(2tiVG), \ e = 0.2. 

In a similar vein, if we define = 2Ek/Eq as the 
ratio of measured gravitational energy to total kinetic 
energy Ek — \ J pv 2 d 3 r, we find for our data at t = 



a E = 0.64. 



(6) 



2.3. 



Physical Scaling 

Although these simulations are scale free, we use a box 
size of 10 pc and density of 300cm -3 for the analysis 
presented here. This gives a resolution on the finest zones 
of 1000 AU and a mean mass density po = 1.2 x 10~ 21 
g. The sound speed is set to be 0.2 kms -1 , giving a 
temperature of ss 10 K. Mean and RMS magnetic fields 
are 0.6 and 2.7 /zG, respectively, with the mean along the 
z direction. These units give a total mass of 1.2 x 10 4 M Q 



Rescaling can be found with the following relations: 



iff 
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i{xT l/2„-l/2 



32Gp 
P 

M oc T 3/2 p- 1/2 
c s oc T 1 ! 2 
Boer 1 / 2 /? 1 / 2 . 
2.4. Gore definition 



(7) 

(8) 
(9) 
(10) 
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We select core s with the algorithm described by 
ISmith et al.l (120091 ) which selects topologically connected 
density isosurfaces at a spacing of Sp — pi/pi—i that have 
no substructure; that is, a contour at Pi-\ that contains 
zero or one cont our above pj. This is a si milar definition 
to tha t used bv lPadoan et al.l (|2007|) and I Schmidt et all 
(|2010[) , both of whom reported no significant effects due 
to variations in dp, provided bp < 1.16. In our simula- 
tion, dp = 1.12. A core is determined to be bound if its 
gravitational energy is at least twice as large as the sum 
of kinetic, thermal, and magnetic energies. All analysis 
has been pe rformed with the AMR analysis package yt 
(lTurkH2008D . 

All analysis in this section takes place at a single snap- 
shot, t = 0.754//. At this time, there are 148 cores that 
match the definition of bound. Figure Q] shows a projec- 
tion of density at this timestep in the left panel, and a 
close up of a portion of the column density map in the 
right panel with core boundaries drawn. All cores drawn 
are associated with topologically isolated, gravitationally 
bound objects, but some have their features washed out 
by projection effects. 



3. DENSITY PDF 



One of the most 
personic turbulence 



prominent 
is the lpj 

of de nsi ties (|Vazauez-Semadenil 
1997allbl: IScalo et al.l 119981: iPassot fc Vazauez-Semadeni 
19981 iNordhmd fc Padoanl 119991 iKlessenl 120001: 
Padoan fc Nordlundl 120021 120091: iFederrath et al.l 



consequences of su- 
- norm al distribution 
119941: iPadoan et al 



2008bD . This has been used to predict several properties 
of star formati on, including the Initial Mass Function 
of sta rs (IMF) (jPadoan fc Nordlundl l200l IPadoan et al. l 
2007 ). brown dwarf frequency (IPadoan fc Nordlundl 
2004 ) and the star formation rate (jKrumholz fc McKed 
20051) . Here we will discuss the PDFs one expects to see 
from isothermal turbulence, and what has been seen in 
our simulations with the inclusion of self gravity. 

The central limit theorem states that the sum of 
a sufficiently large number of uncorrelated events will 
form a Normal, or Gaussian, distribution. A corol- 
lary of this is a sufficiently large number of random 
multiplicative events will form a lognormal distribu- 
tion. This distribution has been experimentally verified 
la rge number of different s i mulations, bo t h pure 



(iVazauez-Semadenil 119941 : IPadoan et al.. 
Kritsuket a l.l 120071 I2 009T) and MHD dLi et al 



hydro 



Lemaster fc Stonell2008l: iKritsuk et aH l2009) 



1997b: 

12004 : 



The log-normal distribution is given by 



P(x)d In x 



1 



V2 



exp 



7T<7 



(lnx — n)' 
2^ 2 



dlnx 



(12) 
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Figure 1. Projection of the full domain at t = 0.75tff (left) and a close up of the squared region, including curves indicating the edges 
of all bound cores in the region (right). 



where x = p/po is the over density, and p = — a 2 / 2 is 
the mean of In a;. For pure hydrodynamical turbulence, 



011(1 + b 2 M 2 ) 



(13) 



2007 



where b has been determined numerically to lie between 
0.3 and 0.4 (IPadoan et alJll997bHFederrath et alll2008bt 

120081 IKritsuk et al l 
The value b has been 
a factor of 3, depend- 
to compressive forcing 



Kritsuk et al 
2010bl iFederrath et al 



Beetz et al 
I2010D 



shown to vary by as much as 
ing on the ratio of solenoid al 

dFederrath et al.H2008bL [20101) 

For driven MHD turbulence. iLemaster fc Stone) (|2008l ) 
find that 



CLM08 



= y/\ - 0.72 In [1 + 0.5.M 2 ] + 0.20| 



(14) 



and is insensitive to the magnetic field strength. 

Figure [2] shows density PDF P(p) for two snapshots: 
at t = (left) and t — 0.75t// (right). In both plots, the 
solid line is the measured PDF and the dashed line is the 
fit to a lognormal. TableQ]shows the fit parameters. The 
addition of self gravity causes the PDF to widen and the 
mean to decrease. We find that b = 0.3 f or our initial con- 
ditions, similar to IKritsuk et al.l (|2007l ). but b = 0.5 for 
the collapsed sna pshot; we also find t hat th e dispersion, 
clmo8, found by ILemaster fc Stond (|2008l ) is in better 
agreement with the collapsed state of the simulation. 

When turbulence simulations are performe d in th e 
prese n ce of sel f grav i ty, several authors (IK esse 
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20001: iSlvz et a .1 120051 : IVazauez-Semadeni et all 
Federrath et al.l l2008at IKritsuk et al.1 l2010a| ) find that 
the log-normal PDF u nderestima t es the high density tail 
of the measured PDF. ISlvz et al.l ([20051 ) fit the h igh den- 
sity t ail to a power law with index of —1.5. iKlessenl 
(|2000l ) does not fit a power law, but the resolution of 
those simulations is much low e r tha n what we present 
here. FVazaucz-Scmadc ni et al.l ()2008l ) mention the exis- 
tence of a power law, but say no thing further. Observa- 
tionally. lKainulainen et all (|2009D found power law wings 
in column density distributions of active star forming re- 
gions, and a similar high density power law has been seen 
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Table 1 

Fit parameters. 

o X 2 



0.75 8.57 



-0.80 
-1.86 



1.35 
1.74 



4.6 x 10" 6 
1.6 x 10" 6 



0.27 
0.52 



°LM08 



1.67 
1.65 



in Aquila (Philippe Andre, private communication.), Our 
w ork is the first reporte d case in an MHD simulation. 

IKritsuk et al.l (|2010a[) find an extended two part power 
law, with index —1.7 at intermediate densities, and —1 
at high densities. They provide the first explanation of 
this power law, associating it with a self-similar singu- 
lar isothermal sphere. In such a sphere, p cx r~ 2 , thus 
V(p) cx p- 3/2 . 

Figure [3] shows both volume-weighted PDF P{p) and 
mass- weighted PDF M(p), with the power law fits of 
— 1.64 and —0.64, respectively, in the range of p = 
10 — 1000. The dashed line in Figure[3]is the same curve 
as the solid line in Figure[2l but here with the mid to high 
density power law emphasized. This power law breaks 
down above a density of 1000, likely due to resolution 
effects, as the maximum resolvable densit y in the simula- 
tion, according to the Truelove condition (|Truelove et al.l 
Il997h . is a p/p Q = 1623. The exponent we find is quite 
close to the —1.7 found bv IKritsuk et al.1 ()2010al ). and in 
reasonable agreement with the —1.5 they predict fr om a 
singular isothermal sphere. iFederrath et al.l (|2008aD also 
measure the density PDF from a driven, self-gravitating 
turbulence simulation (without magnetic fields) with a 
set of Lagrangian tracer particles. The tracer particles 
in their work also follow a power law tail, with an index 
of —0.6 ± 0.1, (Fedderath 2010, private communication) 
consistent with our mass- weighted M(p). 

4. MAGNETIC PDF 

Figure H] shows the volume- weighted PDF of the mag- 
netic field at t = in a sem i-log plot. In agreement with 
IPadoan fc Nordlundl (Il999l) . we find a roughly exponen- 
tial tail at high field strength. This intermittent distri- 
bution in the magnetic field is caused by the large field 
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Figure 2. Density PDFs of the initial conditions (left) and after t = 0.75tff (right). Both are fit to lognormal distributions (dashed 
lines). 
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Figure 3. Volume (left) and mass (right) weighted density PDFs at i = 0.75tff, with power law fits. 



amplifications by the strong, three dimensional compres- 
sions made possible by the large kinetic energy, relative 
to the magnetic energy, of the super- Alfvenic turbulence. 

The left panel of Figure [5] shows the probability den- 
sity function for the magnetic field strength, P(B) for 
t = (solid line) and t — 0.75t// (dashed line), here in 
a log- log plot. The most significant aspect of this fig- 
ure is the strong power law tail after the collapse has 
evolved. Gravitational collapse amplifies the peak mag- 
netic field strength by 2 orders of magnitude (in fact a 
factor of 320), and creates a prominent power law tail. 
This powerlaw tail is fit by P{B) cx B x , with x = —2.74. 



If one naively takes P(p) 



cx p 



-1.5 



as expected from a 



singular isothermal sphc: 
(as in other simulations 



re (as in Secti on |3|) and p oc B 2 



(IBertoldi fc McKedll992 



_i et al.ll2004[) a nd observations 
Crutcherlll999l) of dense cores) , 



one arrives at x = —3, which is quite close to the value 
we find here. Details will be discussed in the next Sec- 
tion, where we measure the relationship between B and 
p in our simulation. 



The right panel of Figure [5] shows the mass- weighted 
PDF, M{B) for the same snapshots as the left panel. 
If we fit a power law to the same field strength range 
that was used for P(B), we find M(B) cx B m , where 
rn = —0.4. However, the power law in M(B) is not 
nearly as well defined as for P{B). This is due two re- 
lated effects: the power law relation between p and B 
is less well defined at field strengths above B = 100/iG 
(see Figure [7] in the next section); and the power law in 
density breaks down above p/po > 1000, likely due to 
resolution effects (see Figure [3]). 

5. FIELD STRENGTH VS. DENSITY 

Figure[6]shows a contour plot of magnetic field strength 
vs density, colored by fraction of mass in each (B, p) 
bin. The left pa nel shows t = 0, b e fore t he action 
of gravity. As in iPadoan fc Nordlundl (|1999l ). the up- 
per envelope is matched by a power law roughly of the 
form B cx p 0A , and the scatter is quite large in both 
B and p. The large scatter in field strength is due to 
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Figure 4. Semi-log plot of P(B) at t = 0. A roughly exponential 
tail can be seen at high field strength, consistent with earlier super- 
Alfvcnic turbulence simulations. 

the fact that only the component of the field perpendic- 
ular to a shock is amplified, but due to the weak na- 
ture of the magnetic field, the relat ive orientations of B 
and th e shock are not correlated. IPadoan fe Nordlundl 
(|1999( ) demonstrated that in models with larger mag- 
netic field strength (sub-Alfvenic turbulence), the flow 
is constrained to move primarily along the field lines, 
significantly decreasing the field strength amplification 
power of the turbulence, and decreasing the scatter in 
the magnetic field strength. The right panel in Figure 
[5] shows the B — p relation at t = 0.75i/ /, after gravity 
has taken effect. Several features are noticeable. Both B 
and p are amplified by almost four orders of magnitude 
by the gravity, and the initial scatter in the distribution 
remains imprinted on the self gravitating distribution. 
The relation gets more shallow at p/p n > 10 3 , which is 
the same density at which the power law slope seen in 
P{p) changes. As before, this is likely due to resolution 
effects. Another feature of figure [6] is that the upper en- 
velope slope increase sharply at p/po > 5. Above this 
density, gravitational collapse has overtaken the dynam- 
ics, as evidenced by the onset of the power law in the 
density PDF. As these motions are three dimensional 
contractions, rather than the one dimensional compres- 
sion due to shocks, the amplification is stronger than 
what is possible from the turbulence alone. 

Figure [7] shows a mass weighted average B as a func- 
tion of density bin, which has been fit to a power law be- 
tween over densities of 10 and 1000. We find B oc p 0A8 . 
This behavior has been predicted or ob s erved by several 
other authors: iFiedler fc Mouschoviasj (|1993[ ) predicted 
B oc p K , with k = 0.44 — 0.5, though this was done 
in a quas i-static collapse model w ith ambipolar diffu- 
sion; both lBertoldi fc McKed (|1992D and lCrutcherl (|1999fl 
found a similar result in dense molecular cores and in- 
terpre ted it as constancy of the Alfven speed; iLi et al.l 
(|2004f) found this behavior for central density and mag- 
netic field strengths. Figure |5] demonstrates that this 
behavior is endemic to the entire collapse process, not 
just the high density collapsed objects. 

6. MASS TO FLUX RATIO 



Figure [5] shows line-of-sight magnetic field strength, 
i?ios, vs column density, N, for three populations 
of cores: CN Zeem an splitting measurements of 
iFalgarone et al.1 (l2008f): OH Zeema n splitting measure- 
ments of iTroland fc Crutcherl ()2008| ): and bound objects 
in our simulation. Color shows A^ = \J Eq / Eb , out- 
adopted proxy for mass to flux ratio. The left plot shows 
t = 0, the right plot shows t = 0.75i//. In the sim- 
ulated points, the line of sight is taken along each of 
the three coordinate axes, so that there are three points 
in the plot for each simulated core; B\ os is the density 
weighted average of each field quantity, and N = M/A, 
the total computed mass of the core divided by the area 
projected along that axis. The observational points are 
currently the best data available to relate magnetic field 
strength to mass, and some of the only magnetic field 
measurements for high density protostellar cores. The 
ratio between these quantities is often used as a proxy 
for mass-to-flux ratio A as 



A 



N(H 2 ) VG 



(•<]> 



(15) 



where c 



(1/2,1/3) is a geometrical correction 



for (spherical, sheet) projection effects, and c$ 
(0.12, l/27r) is a correction found numerically for equi- 
librium configuratio ns of a (sphere, sheet), respectively 
(jBourke et al.ll200TI ) . Here, we directly compare B\ os and 
N(H2). This allows us to eliminate the need for cither 
correction factor and compare the results of our model 
directly to the observations. 

We find that our cores and the observed cores have 
almost the same distribution in the B\ oa -N space. This 
shows that the early evolution of prestellar cores is well 
reproduced with isothermal super Alfvenic turbulence 

and self-gravit y. Simulations without self gravity by 

iLunttila et al.l (|2008l) successfully reproduce the lower 
density OH measurements, but lack the density range 
to reproduce the CN observations. This is also seen in 
the left panel of Figure [8l which shows the initial time 
t = 0, at which only the effects of turbulence are felt by 
the gas. The inclusion of self gravity and the large range 
of density scales allowed by AMR reproduces the higher 
density CN points. 

The color of the simulated points corresponds to Xe = 
J 'Eg I 'En, non -spherical analog of A = (M/$)/(M/$) c . 
ILi et al.l (|2004[ ) found that A > 10 for all cores in ques- 
tion. We find that Xe > 1 for all objects, which, using 
the spherical guide, is analogous to A > 5, so 

our results are in reasonable agreement with theirs. 

The simulated points at t = 0.75i// in figure|8]are best 
fit by a power law 



Bins =N° 



57 



(16) 

Collapse that preserves mass-to-flux would have an ex- 
ponent of 1, by equation [15j Since ideal MHD preserves 
M/$ along a flux tube, this indicates that flow along the 
field lines must be r esponsible for some o f the dynam- 
ics. As discussed by IPadoan fc Nord lund (1999|), this is 
due predominantly to flow along the magnetic field lines. 
This flow is due predominantly to kinematic alignment 
between the velocity and magnetic field, wherein the ki- 
netic energy stretches the magnetic field, hence aligning 
the two. The relationship in equation [16] is also expected 
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Figure 5. Magnetic field PDF P(\B\) (left) and mass-weighted PDF M(\B\) (right) for t = (solid line) and t = 0.75t ff (dashed line) 
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Figure 6. Magnetic field strength vs. density for all zones in the simulation for t 
total mass fraction in each (B, p) bin. 



(left) and t = 0.75i// (right). Color field shows 



from figure demonstrating that the local properties of 
the core are dictated by the global flow properties. 



7. CORE MASS FUNCTION 



One of the open questions in star formation i s the ori- 
in of the stellar initial mass function (IMF). ISalpeterl 
19551 ) first measured this and fit it to a power law, 



o 



0.03(^)W. 
= -2.35 



(17) 
(18) 



This fit was done between 1 and 10 M©. The exact 
value of the expone nt in the 1 to 10 M range is still 
under investigation (|Scaloll2005h . though recent measure- 
ments give the range of a to be between -2.3 and -2.8. 

It has been proposed that the IMF and CMF are di- 
rectly related to one another, either directly (jMotte et al.l 



119981 ) or with some fraction of each core lost in t he final 
collapse and accretion phase (|Enoch et al.l l2008h . This 
implies that the IMF is determined by the global or large 
scale processes of star formation, in our model the com- 
bi ned effects of turbulence an d gravity, as in the model 
of iPadoan fc Nordlundl ()2002f ). Alternative models have 
the IMF set by local physics, once protostars have formed 
within the prestellar condensations. Th ese models in- 
clude the competitive accretion model of IBonnell et al.l 
(|2001[) , wherein the population of neighboring protostars 
in fluences the final m ass of any giv en star; and models 
of iShu et all (|1987t) or iMversI (|2010[) . where protostellar 
outflows halt or slow the inflow of gas onto the protostar. 

Figure [S] shows the mass distribution for all bound 
cores. The fit to the high end of the distribution is 
n(M) cx M- 2A±0 - 6 . The fit was performed by fitting 
a power law between the peak and the highest bin for 
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Figure 7. Mass weighted average B vs. p (dotted line) and fit to 
the range where the density PDF is also a clear power law (solid 
line). The fit is a power law with exponent 0.48. 

a succession of bins between 5 and 25. We find good 
agreement between our slope and the IMF slopes men- 
tioned ab ove, and the s l ope of — 2.3±0.6 measured for the 
CMF bv lEnoch et al.1 (|2008l ). Further agreement with 
the CMF is seen in figure [TU1 which shows the cumu- 
lative mass function iV(> M) = J M n(M)dM for our 
data (black line) and the prestellar cores from Perseus , 
Ophiuchus and Serpens presented in lEnoch et al.l (|2009f ) 
(grey line). Here we have scaled both populations to the 
Bonnor Ebert mass, 



M BE = 



l-18c 4 
G 3 / 2 Po /2 ' 



The observed points used a Bonnor Ebcrt mass of 
1.5M Q , which corresponds to a background density of 
w 9000cm -3 at 10 K. This is somewhat higher than the 
mean density in these clouds, but not unreasonable for 
the ambient density immediately surrounding the cores, 
which have mean densities of « 10 5 cm -2 . The simulated 
cores used the mean density in the box, which in the 
scaling used throughout this paper gives M BE = 10Af Q , 
though strictly speaking this is a free value. The observa- 
tional data are all prestellar cores from the 1.1 mm Bolo- 
ca m survey of Perseus , Serpens and Ophiuchus presented 
in lEnoch et al.l (|2006t) . lYoung et al.l (|2006h . lEnoch et al.l 
(|2007|) . respectively, that do not have an associated in- 
frared source in the Spitzer c2d catalog. The majority 
of these objects can reasonably be assumed to be self 
gravitating, based on comparisons of the Perseus cores 
to kinematic informati on from the molecular li ne survey 
of the same region by iRosolowskv et ail ([2008D . As our 
simulation only attempts to model the prestellar core 
phase of star formation, and not the formation of the 
actual star itself, this sample of objects is the best ob- 
servational counterpart for comparison. 

8. CONCLUSIONS 

In this work, we present density and magnetic field 
distributions for a super Alfvenic turbulence simulation 
with self gravity. The simulation was run with the AMR 
extension of Enzo described by iCollins et al.l (|2010l ). 



allowing us unprecedented spatial resolution. Super- 
Alfvcnic turbulence has been proposed as the primary 
mechanism for star formation, pro viding good explana- 
tions of the star formation rat e (jKrumholz fe McKed 
120051: IPadoan fc Nordlundl 120091) and initial mass func- 
tion (jPadoan fc Nordlundl 12002) . Here we provide two 
checks of this model against observations, and explore 
deviations from the predictions of super- Alfvenic turbu- 
lence caused by the addition of self gravity. 

We find in Sections[3]and[3]that power law tails develop 
for both high density and high magnetic field in volume 
and mass weighted PDFs, P{p) oc p~ 167 and P{B) oc 
B~ 2 74 , respectively. The volume-weighted density PDF 
is consistent with the prediction of a singular isothermal 
sphere (SIS) (jKritsuk et al.H2010aD P(psis) oc p~ 15 . 

The relationship between the magnetic field and the 
density also shows a power law behavior B oc p 0AS 
throughout th e gas, consistent with the findings of 
iLi et al.l (|2004f ) , who found a similar behavior in the peak 
density/field relation in cores. This then allows us to ex- 
plain the magnetic PDF, by combining this result with 
the density PDF. 

Gravitationally bound cores found in our simula- 
tion were compared against several observational sur- 
veys. Comparisons wit h the most recent Z e eman split- 
ting measurements of iTroland fc Crutcherl (|2008| ) and 
iFalgarone et al.l (|2008f ) show that the mass-to-flux ra- 
tio in our simulations agrees in value and behavior with 
those found observationally. The relationship between 
field strength and column density is fit to a power law, 
B oc N - 57 , demonstrating that significant mass-to-flux, 
thus magnetic support, is lost due to motion along the 
field lines. 

Comparing our core mass function (CMF) to that of 
prestellar cores in lEnoch et al.l (|2008l) . we again find ex- 
cellent agreement. A slope of 2.1 ± 0.6 agrees with their 
fit value of 2.3 ± 0.6, and cumulative mass distributions 
line up almost identically. The relatively good match 
between the observed CMF and the observed IMF indi- 
cates that the IMF is determined well before the onset 
of nuclear burning, at a relatively low (compared to the 
protostar) density. A mul tiplicative off s et of > 1/4 is 
seen between the cores of lEnoch et all {2008) and the 
observed IMF, indicating that as much as 3/4 of the 
mass is lost in the final collapse phase. However, this 
is only an upper limit to the lost fraction, as the peak 
of the observed CMF is heavily influenced by its com- 
pleteness limit. Bound cores in our simulation agree with 
the observed CMF extremely well, indicating that super- 
Alfvenic turbulence and gravity are primarily responsible 
for the structure of the mass distribution of the CMF, 
and ultimately the IMF. 

The authors would like to thank A. Kritsuk for his 
input and many useful discussions. The authors would 
like to acknowledge financial support from NSF grants 
AST0808184 and AST0908740, and computational re- 
sources provided by the National Institute for Computa- 
tional Sciences under LRAC allocation MCA98N020 and 
TRAC allocation TG-AST090110. 
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